Play building blocks on population distribution of multilevel superconducting flux 

qubit with quantum interference 
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Recent experiments on Landau-Zener interference in multilevel superconducting flux qubits re- 
vealed various interesting characteristics, which have been studied theoretically in our recent 
work[PRB 79,094529, (2009)] by simply using rate equation method. In this note we extend this 
method to the same system but with larger driving amplitude and higher driving frequency. The 
results show various anomalous characteristics, some of which have been observed in a recent ex- 
periment. 

PACS numbers: 74.50.-|-r, 85.25.Cp 



In this note, we show the anomalous characteristics 
in population distribution of superconducting flux qubit 
under driving fields with large amplitude(A) and high 
frequency(a;). For a field-driving multilevel flux qubit 
with diabatic quantum states |i,i)(left well, i=0,l,2...) 
and |j, i?)(right well, j=0,l,2,...),— the Landau-Zener 
transition rate between \i,L) and \ j,R) is>^ 
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where Ay is the avoided crossing between states \i,L) 
and \j, R), ,eij is the dc energy detuning from the corre- 
sponding avoided crossing Ay , r2 = 1 /T2 is the dephas- 
ing rate and Jn{x) are Bessel functions of the first kind 
with the argument x — A/lo. 

The time evolution of population for state \i,L) can be 
described bj*^ 
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where Ti^ii is the intrawell relaxation rate from \i,L) 
to \i' , L) , and Ti^j is the interwell relaxation rate from 
\i, L) to I J, R). The time evolution of population for state 
\j, R) can be obtained in the same way. In the station- 
ary case, we have Pi^L{R) = 0. The qubit population 
distribution can be easily obtained by 



Based on Eq.([T][3]), we can play building blocks. Eq.((T]) 
serves as the basic building block, Eq.([2]) serves as the 
rule of the game, and Eq.Q shows the total pattern. 
By adjusting the frequency w and dephasing rate r2 
in Eq.([T]), we can choose different kinds of building 
blocks. Therefore, kinds of interesting patterns can be 
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Field Amplitude 

FIG. 1: Calculated qubit population versus flux detuning and 
driving amplitude. The parameters are from experiments.— 
The dephasing rate is r2/27r = 0.5GHz and the driving fre- 
quency is u!/2tt — 160MHz. 



constructed. As shown in Fig. 1, the diamond pat- 
terns observed in recent experimentsi>^ are well obtained 
based on Eq.([l][3]) with experimental parameters. By in- 
creasing the driving frequency continually, we can obtain 
more interesting patterns. As shown in Fig. 2, with 
the parameters in a recent experiment where the driving 
frequency is much higher(a;/27r~10GHz), we can obtain 
various anomalous patterns, some of which have been 
demonstrated by experiments.— 



^ D. M. Berns, M. S. Rudner, S. O. Valenzuela, K. K. 

Berggren, W. D. Oliver, L. S. Levitov, and T. P. Orlando, 

Nature (London) 455, 51 (2008). 
^ W. D. Oliver and S. O. Valenzuela, Quantum Inf. Process. 

8, 261 (2009). 



^ X. Wen and Y. Yu Phys. Rev. B 79, 094529 (2009). 
* Y. Wang, S. Cong, X. Wen, C. Pan, C. Sun, J. Chen, L. 
Kang, W. Xu, Y. Yu, and P. Wu, arXiv:0912.0995. 




static Flux Detuning ( m^o) 



FIG. 2; Calculated qubit population versus flux detuning 
and driving amplitude. The simulation parameters are from 
experiments.* The driving frequencies(tj/27r) are (a)5GHz, 
(b)8GHz, (c)llGHz, (d)13GHz, (e)15GHz, and (f)17GHz, re- 
spectively. 



